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MATCHED VIEW WEIGHTING METHODS AND 
APPARATUS USING MULTIPLE TILTED 
RECONSTRUCTION PLANES 

BACKGROUND OF THE INVENTION 

[0001] This invention relates generally to computed tomography 
(CT) and more particularly to matched view weighting methods for computed 
tomography scanners using multiple tilted reconstruction planes. 

[0002] A single tilted-plane-based reconstruction (sTPR) algorithm 
has been proposed to combat cone beam and helical artifacts in helical multi-detector- 
row CT. Each tilted reconstruction plane (RP) in sTPR uses two freedoms of rotation. 
A virtual x-ray source trajectory (VT) is defined as the intersection of a RP with the 
cylinder containing an actual helical x-ray source trajectory (AT) to reconstruct an 
image on each tilted RP. The error between VT and AT (namely trajectory error) is 
globally minimized and symmetric if the central intersection point between a RP and 
an AT is taken as the origin. However, noise characteristics of sTPR may be 
compromised because only short scan (180° + fan angle) projection data are used to 
reconstruct an image on a tilted RP. To improve the noise characteristics of sTPR, a 
multiple-tilted-plane-based reconstruction (mTPR) algorithm has been proposed, in 
which each tilted RP uses three freedoms of rotation. As a result, the trajectory error 
in mTPR is no longer globally minimized and symmetric. 

BRIEF DESCRIPTION OF THE INVENTION 

[0003] In one aspect, a method for Computed Tomography scanners 
using multiple tilted reconstruction planes is provided. The method includes view 
weighting data in accordance with 
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is a central view angle, y is a fan angle, y m is a maximum half fan angle of an x-ray 
beam, and w b and w t are each weights no greater than 1. 

[0004] In another aspect, a Computed Tomography (CT) system 
includes an x-ray source, a detector positioned to receive x-rays emitted from the 
source, and a computer operationally coupled to the source and the detector. The 
computer is configured to receive data from a scan of an object, and weight the 
received data in accordance with 

w b 0<P<2y m -2y 
Z(/3,/) = \ 1 2y m -2y<p<x-2y <o(P,Y) = 3 X \P,r)-2x\P,y) where 0 

w, n -2y < P <7t + 2y m 

is a central view angle, y is a fan angle, y m is a maximum half fan angle of an x-ray 
beam, and w b and w t are each weights no greater than 1 . 

[0005] In another aspect, a method of obtaining and using weights is 

provided. The method includes setting w b = — and w t = — with 

e,+e b e,+e b 

e b = j] AT(fi) - VT(j5) | dp and e, = J| AT(p) - VT(0) | dfi wherein AT(/8) is an 

0 * 

actual helical x-ray source trajectory (AT) of 13 and VT(/3) is a virtual x-ray source 
trajectory (VT) of j3, and weighting image data with w b and w t prior to filtering and 
backprojecting the image data to reconstruct an image. 

[0006] In still another aspect, a computer readable medium encoded 
with a program configured to instruct a computer is provided. The program is 
configured to instruct the computer to receive data from a scan of an object, and 
weight the received data in accordance with 

w b 0<p<2y m -2 r 
Z(P y y) = \ 1 2y m -2y<p<n-2y o(P,y) = 3 X 2 (P,y) -2 X 3 (P,y) where p 
w, K-2y < p <k + 2y m 
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is a central view angle, y is a fan angle, y m is a maximum half fan angle of an x-ray 
beam, and w b and w t are each weights no greater than 1 . 

BRIEF DESCRIPTION OF THE DRAWINGS 

[0007] Figure 1 is a pictorial view of a CT imaging system 

embodiment. 

[0008] Figure 2 is a block schematic diagram of the system 
illustrated in Figure 1 . 

[0009] Figure 3 illustrates that a trajectory error in sTPR is globally 
minimized and symmetric over view angle. 

[0010] Figure 4 illustrates a 2D intensity distribution of a view 
weighting function. 

[0011] Figure 5 illustrates a weighted trajectory error (mm) by 
applying the view weighting function shown in Figure 4 on the trajectory error shown 
in Figure 3. 

[0012] Figure 6 illustrates a trajectory error in mTPR increases and 
no longer symmetric about the central attachment point. 

[001 3] Figure 7 graphically illustrates a view weighting function. 

[0014] Figure 8 illustrates a weighted trajectory error (mm) by 
applying the view weighting function shown in Figure 7. 

[0015] Figure 9 graphically illustrates a view weighting function. 

[0016] Figure 10 is a three dimensional perspective view of a helical 
body phantom (HBP). 

[0017] Figure 11 illustrates two tomographic image of the HBP 
phantom shown in Figure 10 reconstructed on a RP, in which the matched and non- 
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matched view weighting functions were applied respectively before filtering and 
backprojection. 

DETAILED DESCRIPTION OF THE INVENTION 

[0018] Herein described are methods and systems to compensate for 
the above described asymmetry in trajectory error in mTPR, utilizing a matched view 
weighting (MVW) (described in detail below), which is optimal from the perspective 
of both suppressing cone beam and helical artifacts and maintaining noise 
characteristics. 

[0019] The apparatus and methods are illustrated with reference to 
the figures wherein similar numbers indicate the same elements in all figures. Such 
figures are intended to be illustrative rather than limiting and are included herewith to 
facilitate explanation of an exemplary embodiment of the apparatus and methods of 
the invention. 

[0020] In some known CT imaging system configurations, an x-ray 
source projects a fan-shaped beam which is collimated to lie within an X-Y plane of a 
Cartesian coordinate system and generally referred to as an "imaging plane". The x- 
ray beam passes through an object being imaged, such as a patient. The beam, after 
being attenuated by the object, impinges upon an array of radiation detectors. The 
intensity of the attenuated radiation beam received at the detector array is dependent 
upon the attenuation of an x-ray beam by the object. Each detector element of the 
array produces a separate electrical signal that is a measurement of the beam intensity 
at the detector location. The intensity measurements from all the detectors are 
acquired separately to produce a transmission profile. 

[0021] In third generation CT systems, the x-ray source and the 
detector array are rotated with a gantry within the imaging plane and around the object 
to be imaged such that the angle at which the x-ray beam intersects the object 
constantly changes. A group of x-ray attenuation measurements, i.e., projection data, 
from the detector array at one gantry angle is referred to as a "view". A "scan" of the 
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object comprises a set of views made at different gantry angles, or view angles, during 
one revolution of the x-ray source and detector. 

[0022] In an axial scan, the projection data is processed to 
reconstruct an image that corresponds to a two dimensional slice taken through the 
object. One method for reconstructing an image from a set of projection data is 
referred to in the art as the filtered back projection technique. This process converts 
the attenuation measurements from a scan into integers called "CT numbers" or 
"Hounsfield units", which are used to control the brightness of a corresponding pixel 
on a cathode ray tube display. 

[0023] To reduce the total scan time, a "helical" scan may be 
performed. To perform a "helical" scan, the patient is moved while the data for the 
prescribed number of slices is acquired. Such a system generates a single helix from a 
fan beam helical scan. The helix mapped out by the fan beam yields projection data 
from which images in each prescribed slice may be reconstructed. To further improve 
the performance of the CT system, multi-slice CT systems are built. In such systems, 
multiple projections are acquired simultaneously with multiple detector rows. For 
example, by extending the x-ray width in the z-axis forming a cone beam, cone beam 
helical scans are performed. 

[0024] Reconstruction algorithms for helical scanning typically use 
helical weighing algorithms that weight the collected data as a function of view angle 
and detector channel index. Specifically, prior to a filtered backprojection process, 
the data is weighted according to a helical weighing factor, which is a function of both 
the gantry angle and detector angle. The weighted data is then processed to generate 
CT numbers and to construct an image that corresponds to a two dimensional slice 
taken through the object. 

[0025] As used herein, an element or step recited in the singular and 
preceded with the word "a" or "an" should be understood as not excluding plural said 
elements or steps, unless such exclusion is explicitly recited. Furthermore, references 
to "one embodiment" of the present invention are not intended to be interpreted as 
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excluding the existence of additional embodiments that also incorporate the recited 
features. 

[0026] Also as used herein, the phrase "reconstructing an image" is 
not intended to exclude embodiments of the present invention in which data 
representing an image is generated but a viewable image is not. Therefore, as used 
herein the term "image" broadly refers to both viewable images and data representing 
a viewable image. However, many embodiments generate (or are configured to 
generate) at least one viewable image. 

[0027] Additionally, although described in a medical setting, it is 
contemplated that the benefits of the invention accrue to all CT systems including 
industrial CT systems such as, for example, but not limited to, a baggage scanning CT 
system typically used in a transportation center such as, for example, but not limited 
to, an airport or a rail station. 

[0028] Referring to Figures 1 and 2, a multi-slice scanning imaging 
system, for example, a Computed Tomography (CT) imaging system 10, is shown as 
including a gantry 12 representative of a "third generation" CT imaging system. 
Gantry 12 has an x-ray source 14 that projects a beam of x-rays 16 toward a detector 
array 18 on the opposite side of gantry 12. Detector array 18 is formed by a plurality 
of detector rows (not shown) including a plurality of detector elements 20 which 
together sense the projected x-rays that pass through an object, such as a medical 
patient 22 between array 18 and source 14. Each detector element 20 produces an 
electrical signal that represents the intensity of an impinging x-ray beam and hence 
can be used to estimate the attenuation of the beam as it passes through object or 
patient 22. During a scan to acquire x-ray projection data, gantry 12 and the 
components mounted therein rotate about a center of rotation 24. Figure 2 shows only 
a single row of detector elements 20 (i.e., a detector row). However, multi-slice 
detector array 1 8 includes a plurality of parallel detector rows of detector elements 20 
such that projection data corresponding to a plurality of quasi-parallel or parallel slices 
can be acquired simultaneously during a scan to enable cone beam helical scans. 
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[0029] Rotation of gantry 12 and the operation of x-ray source 14 are 
governed by a control mechanism 26 of CT system 10. Control mechanism 26 
includes an x-ray controller 28 that provides power and timing signals to x-ray source 
14 and a gantry motor controller 30 that controls the rotational speed and position of 
gantry 12. A data acquisition system (DAS) 32 in control mechanism 26 samples 
analog data from detector elements 20 and converts the data to digital signals for 
subsequent processing. An image reconstructor 34 receives sampled and digitized x- 
ray data from DAS 32 and performs high-speed image reconstruction. The 
reconstructed image is applied as an input to a computer 36 which stores the image in 
a storage device 38. 

[0030] Computer 36 also receives commands and scanning 
parameters from an operator via console 40 that has a keyboard. An associated 
cathode ray tube display 42 allows the operator to observe the reconstructed image 
and other data from computer 36. The operator supplied commands and parameters 
are used by computer 36 to provide control signals and information to DAS 32, x-ray 
controller 28 and gantry motor controller 30. In addition, computer 36 operates a 
table motor controller 44 which controls a motorized table 46 to position patient 22 in 
gantry 12. Particularly, table 46 moves portions of patient 22 through gantry opening 
48. 

[0031] In one embodiment, computer 36 includes a device 50, for 
example, a floppy disk drive, CD-ROM drive, DVD drive, magnetic optical disk 
(MOD) device, or any other digital device including a network connecting device such 
as an Ethernet device for reading instructions and/or data from a computer-readable 
medium 52, such as a floppy disk, a CD-ROM, a DVD or an other digital source such 
as a network or the Internet, as well as yet to be developed digital means. In another 
embodiment, computer 36 executes instructions stored in firmware (not shown). 
Computer 36 is programmed to perform functions described herein, and as used 
herein, the term computer is not limited to just those integrated circuits referred to in 
the art as computers, but broadly refers to computers, processors, microcontrollers, 
microcomputers, programmable logic controllers, application specific integrated 
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circuits, and other programmable circuits, and these terms are used interchangeably 
herein. 

[0032] Figure 3 illustrates that the trajectory error in sTPR is globally 
minimized and symmetric over view angle. More specifically, Figure 3 illustrates a 
variation of trajectory error (mm) as a function over central view angle P (°) in sTPR 
(detector row number: 64, detector row width: 0.625 mm, and normalized pitch of 
1 .0). A few view weighting functions, most of which are symmetric about a central 
view angle, have been proposed to reconstruct short scan data. For instance, 

P 



°' 5 y m -y 0<P<2 Ym -2y 

1 ' 2y m -2y<j3<7r-2y (1) 

0.5 * + 7m — — n - 2y < p < n + 2y m 



Ym +7 

co(P,y) = 3 Z 2 (P,r)-2z\P,r) (2) 

where (5 is the central view angle, y the fan angle, and y m the maximum half fan angle 

of an x-ray beam (see Figure 4 which illustrates a 2D intensity distribution of the view 
weighting function given in (1) and (2) for fan beam short scan data wherein white 
pixels correspond to 1 .0, gray to 0.0-1 .0 and black to 0.0). This view weighing 
function appears apt for sTPR because its symmetric distribution matches the 
trajectory error. Figure 5 illustrates a weighted trajectory error (mm) by applying the 
view weighting function shown in Figure 4 on the trajectory error shown in Figure 3 
(detector row number: 64, detector row width: 0.625 mm, normalized pitch 1.0). 

[0033] However, and now referring to mTPR, each tilted RP in 
mTPR makes use of three freedoms of rotation along x, y and z axes respectively. 
The rotations along y and z axes are the same as that in sTPR to fit a helical source 
trajectory, and that along x axis is to improve noise characteristics or dose efficiency. 
Intuitively, due to the rotation along the x axis, each tilted RP in sTPR splits into 
several titled RPs in mTPR, and these RPs deviate from the optimized pose 
determined in sTPR. Consequently, the trajectory error in mTPR increases and is no 
longer symmetric about the central attachment point as shown in Figure 6, where 
Figure 6a illustrates a distribution corresponding to the RP rotating along x axis 
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toward -z; and Figure 6b illustrates a distribution corresponding to the RP rotating 
along x axis toward +z; (detector row number 64; detector row width 0.625 mm; 
normalized helical pitch 1 .0; only double tilted reconstruction planes are considered). 
Eventually, the severity of cone beam and helical artifacts are determined by weighted 
trajectory error. A way to compensate for the asymmetry in trajectory error is to use 
an opposite asymmetry in view weighting function. Extreme examples of such 
opposite asymmetric view weighting functions are given by 



,„ N [1 0<B<]t-2y 
wXP,y) = \ (3) 
,KH |p 7r-2y<p<jz + 2y m W 

,„ x fO 0<p<2y m -2y 

n(Ar)= , o o <*< 1 (4) 

[1 2^ m -2y</3<7r + y m 
[0034] Figure 7 graphically illustrates the view weighting function 
given by (3) and (4). Figure 8 illustrates that the view weighting function given by (3) 
and (4) is the best in terms of weighted trajectory error in mTPR, but the worst from 
the perspective of noise characteristics, since all the weighting factors equal to either 
one or zero. 



An optimized view weighting function of mTPR is 



x(P,r) = { 1 



0<P<2y m -2y 
2y m -2y<p<n-2y 
x -2y < p < x + 2y m 



(5) 



a>(P,y) = 3x 2 (P,r)-2% 3 (P,y) 



(6) 



which matches the asymmetric distribution of trajectory error in mTPR. The 
weighting factors w b and w t corresponding to the bottom and upper redundant areas 
respectively are adaptively determined by 



(7) 



e,+e b 



(8) 
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where 

e b = \\AT{P)-VT{P)\dp (9) 

0 

e, = ]\AT(j3)-VT(j3)\dj3 (10) 

[0035] Obviously, there exist discontinuities in the view weighting 
function given in (5) - (10). To reduce or eliminate streak artifacts in reconstructed 
tomographic images caused by those discontinuities, feathering, such as that described 
by C. R. Crawford and K. F. King, "Computed tomography scanning with 
simultaneous patient translation," Med. Phys. 17(6), pp967-982, 1990., may be 
applied while view weighting. Figure 9 graphically illustrates the view weighting 
function given in (5) - (10). 

[0036] Figure 10 is a three dimensional perspective view of a helical 
body phantom (HBP) 60 which was employed to evaluate the imaging performance of 
the herein described matched view weighting function in mTPR. HBP 60 is 
configured to simulate the human thorax under helical scanning. HBP 60 includes a 
plurality of small cylindrical rods which are oblique relative to a gantry's rotation axis 
to simulate the ribs in a human thorax, which would highlight helical artifacts due to 
projection data inconsistence in helical scanning. The x-ray technique used for 
scanning HBP phantom 60 was 120 kVp, 200 mA and 1 .0 sec. per gantry rotation, and 
the scanner possessed 16 detector rows at width 1.25 mm running at normalized pitch 
15/16. A quantitative evaluation of the helical artifacts and noise characteristics was 
conducted by measuring a global helical artifact index (ART + ) as discussed by T. L. 
Toth and P. Simoni, ("A quantitative measure of CT helical Artifacts," Radiology, 
221(P), pp218, 2001.) and the standard deviation (SD + ) within a square containing 20 
by 20 pixels at the isocenter. For simplicity, only double tilted reconstruction planes 
were used at each central view position. The tomographic image on a RP was 
obtained through fan beam, rather than parallel beam, reconstruction algorithm. Both 
w b and w t were varied in step 0.1 respectively in the quantitative evaluation as seen in 
table 1. 
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[0037] Table 1 





w„ 


Sd+ 


Art+ 


1.000 


0.000 


42.214 


14.647 


0.900 


0.100 


41.537 


14.182 


0.800 


0.200 


41.205 


13.525 


0.700 


0.300 


40.942 


13.006 


0.600 


0.400 


40.776 


12.495 


0.500 


0.500 


40.620 


12.278 


0.400 


0.600 


40.661 


11.974 










0.300 


0.700 


40.652 


11.839 


0.200 


0.800 


41.006 


11.948 


0.100 


0.900 


41.104 


12.021 


0.000 


1.000 


41.383 


11.860 


CK 


41.161 


12.538 



[0038] Table 1 illustrates noise characteristics and artifact indexes 
measured in images reconstructed by the mTPR in which the matched view weighting 
function is applied, where CK represents the view weighting function given in Eqs. 
(1) - (2), and the superscript + means tilted RP's 3rd rotation is toward + z axis. The 
artifact index is measured on tilted tomographic images before z-interpolation used to 
obtain tomographic images orthogonal to the gantry rotation axis. It is observed that, 
Art + reaches its maximum at Wb = 0.0 due to a complete mismatch between the source 
trajectory error distribution and the view weighting function. The Art + is inversely 
proportional to the increment in w b and reaches the minimum at w b = 0.7 (w t = 0.3). 
Notice that, as shown in the shadowed row in Table 1, w b = 0.7 (w t = 0.3) is very 
close to the matched view weighting factor w b = 0.608 (w t = 0.392) predicted by Eqs. 
(7) - (10), which is consistent with the expectation on matched view weighting. 
Meanwhile, the artifact index corresponding to the view weighting function given in 
Eqs. (1) - (2) is also listed in Table 1 as a reference (CK), which is significantly larger 
than its counterpart corresponding to the matched view weighting at w b = 0.608 (w t = 
0.392). 
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[0039] The noise characteristics (Sd + ) of the images reconstructed by 
the mTPR on the tilted RP are also listed in Table. 1, where the superscript + means 
tilted RP's 3rd rotation is toward + z axis as before. The standard deviation is 
measured at the isocenter of a tilted tomographic image before z-interpolation. It is 
observed that Sd + reaches the minimum while Wb = w t = 0.5, because the weighting 
factor equaling to 0.5 is the best from the perspective of noise characteristics. 
Moreover, the Sd + is quite stable between w b = 0.5 (w t = 0.5) and w b = 0.7 (w t = 0.3), 
and w b = 0.608 (w, = 0.392) falls into this interval. This is again consistent with the 
expectation on matched view weighting. Meanwhile, the noise characteristics 
corresponding to the view weighting function given in Eqs. (l)-(2) is also listed in 
Table 1 (CK), which is significantly larger than that corresponding to the matched 
view weighting at w b = 0.608 (w t = 0.392). 

[0040] In order to suppress helical artifact as much as possible, only 
the projection data corresponding to an angle range of (180° + fan angle) are used to 
reconstruct an image on tilted reconstruction plane in single/multiple tilted-plane- 
based reconstruction algorithm. The view weighting functions usually used in non- 
tilted-plane-based short scan reconstruction can be used in sTPR because the x-ray 
source trajectory error distribution is symmetric about the central attachment point. 
However, the tilted reconstruction planes in mTPR deviates from the globally 
optimized pose and the resultant x-ray source trajectory error is no longer symmetric 
about the central attachment point. Accordingly, as illustrated above, opposite 
asymmetry in view weighting function is useful to maintain mTPR's capability of 
suppressing helical artifacts while keeping the noise characteristics good. 

[0041] Figure 11 illustrates two tomographic image of the HBP 
phantom reconstructed on a RP by mTPR, in which different view weighting 
functions were applied before filtering and backprojection. Figure 11a illustrates the 
non-matched view weighting given in Eqs. (1) - (2), and Figure lib illustrates the 
matched view weighting given in Eqs. (5) - (6) at wt = 0.392, wb = 0.608. 

[0042] It is observed that the matched view weighting function can 
improve both the noise characteristics and capability of suppression helical artifact 



-12- 



131201 



Express Mail No. EV 339984374 U|S 



significantly. Moreover, as elucidated above, the visual observation in Figure 1 1 of 
improved noise characteristics and helical artifact suppression is reinforced by the 
quantitative results demonstrated in Table 1 . 

[0043] While the invention has been described in terms of various 
specific embodiments, those skilled in the art will recognize that the invention can be 
practiced with modification within the spirit and scope of the claims. 
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